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Abstract. An exact nonsingular solitary wave solution of the Schafer-Wayne 
short pulse equation is derived from the breather solution of the sine-Gordon 
equation by means of a transformation between these two integrable equations. 



The short pulse equation (SPE), which has the form 

U x* = «+H" 3 ) M C 1 ) 

up to scale transformations of its variables, was derived by Schafer and Wayne 
as a model equation describing the propagation of ultra-short light pulses in silica 
optical fibers. In contrast to the nonlinear Schrodinger equation (NLSE) which 
models slowly varying wave trains, the SPE approximates Maxwell's equations in 
the case when the pulse spectrum is not narrowly localized around the carrier fre- 
quency, that is when the pulse is as short as a few cycles of the central frequency. 
Numerical simulations made in % showed that the accuracy of the SPE approxi- 
mation to Maxwell's equations steadily increases as the pulse shortens, whereas the 
NLSE approximation fails to be accurate for ultra-short pulses. 

Like the NLSE, the SPE is an integrable nonlinear partial differential equa- 
tion. Actually, the equation QJ appeared first a long time ago, as one of Rabelo's 
equations possessing a zero-curvature representation with a parameter [3] (we are 
grateful to Prof. E.G. Reyes for this essential reference; see also 0]). Recently, the 
Lax pair of the SPE, of the Wadati-Konno-Ichikawa (WKI) type, was rediscovered 
in [H] , where also the second-order recursion operator of the SPE was presented and 
the chain of transformations 

v{x,t) = (u x + l) ; 
x = w(y, t), v(x, t) = w y (y, t); (2) 
z(y,t) = arccoswy 

was constructed which relates the SPE Q with the sine-Gordon equation (SGE) 

z yt = sinz. (3) 

Very recently, the bi-Hamiltonian structure of the SPE was discovered and the 
corresponding hierarchy and conservation laws were studied in [HI Ej . 

The soliton solutions of the NLSE have played an important role in the recent 
development of fiber-optic communications. For this reason, one would definitely 
like to have exact mathematical expressions for the ultra-short light pulses governed 
by the SPE. In the present paper, we derive exact solitary wave solutions of the 
SPE JJJ from the well-known ones of the SGE J3J by means of the transformation 
£D between these two integrable equations. 
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Schafer and Wayne pQ proved the nonexistence theorem that the SPE does not 
possess any solution representing a smooth localized pulse moving with constant 
shape and speed. Actually, it is easy to find the "wave of translation" solution of 
the SPE directly, in order to see that this solution is singular. Substituting 

u = tt(£), £ = x + ct, c = constant (4) 

into (0 i we get the second-order ordinary differential equation 

(u 2 - 2c) + 2u (u\ + 1) = 0. (5) 

Then, using the integrating factor 2 (u 2 — 2c) Ug, we get from JSJ the first-order 
equation 

(u 2 - 2c) 2 (u 2 + 1) = d, (6) 

where we must set d — Ac 2 for the constant of integration d, in order that the 
asymptotic conditions \u\ — ► and \u^\ — » at |£| — > oo be satisfied. The equation 
© in its form 



uy 4c — u 2 



clearly shows that the sought solution u inevitably reaches the value u = v2c or 
u = — \[2c, where the solution has the singularity — > ±oo. Further integration 
of the equation 10 is possible, and it would bring us to the loop soliton solution of 
the SPE 0). We will obtain this solution soon, but in a different way. 

For simplicity of the subsequent practical use, we can rewrite the transformation 
|2"| in the following equivalent form: 

u(x,t) = z t (y,t), (8) 

{Wy = COS z, 
1 2 ( 9 ) 
w t = - 5 Zf 

The system of two equations for w(y, t) inside I© is compatible if and only if z(y, t) 
satisfies the SGE For any solution z of the SGE 10, the relation JHJ uniquely 
determines it as a function of y and t, and the relation determines x as a function 
of y and t up to an additive constant of integration. This yields a solution u(x,t) 
of the SPE ijljl. given in a parametric form, with y being the parameter. Thus, 
using the transformation |JS}— (JSJ, we can easily derive exact solutions of the SPE 
from the known ones of the SGE. Note that the invariance of the SGE J2J under 
the Lorentz transformation 

y h- > ay, 1 1— > a~ 1 t, z z (10) 

corresponds via (JBJ) © to the invariance of the SPE under the scale transfor- 
mation 

x i — ► ax, t i ► a _1 i, u h-> au, (11) 

where a = constant 7^ 0. Therefore we can put the source solution of the SGE into 
a simpler form by i|l(J|) . thus simplifying the symbolic integration required in JJjJ, 
and then use l|ll|) to generalize the target solution of the SPE. Also we can simplify 
the source solution of the SGE by shifts of y and t, y 1— > y + yo and t t + t : 
a shift of £ in z{y, t) causes the same shift of t in u(x, t), while a shift of y has no 
effect on the target solution of the SPE. 
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Figure 1. The loop soliton solution u(x,t) JT3J) of the SPE (QJ: 

(A) t = -2; (B) t = 6. 



Now let us see how the transformation ©-© works. It is natural to begin with 
the kink solution of the SGE © 8 , in the form 



4 arctan(exp(y + 1)) 



(12) 



simplified by a Lorentz transformation and a shift of y. In the case of z(y, t) given 
by (|12[l. the relations JHJ and © immediately yield 



u = 2 sech(y + t), x — y — 2 tanh(y + t), 



(13) 



where we have fixed the constant of integration in x so that x\ y= t=o = 0. The 
obtained solution u(x,t) l|13|) of the SPE is the loop soliton which moves from 
the right to the left with constant shape and unit speed, as shown in Figure^ This 
multi-valued solution satisfies the equation (JZJ) with c = 1 and has two singularities 
with u x = Ut — » ±oo at the level u = >/2. (Note that, if we chose the antikink 
solution of the SGE as a source solution for the transformation (JSJl- @, the target 
solution of the SPE would be the antiloop soliton which differs from the loop soliton 
in the sign of u only.) We can generalize the solution i|13|) by the transformation 
(|llfl . thus obtaining either a bigger and faster loop or a smaller and slower one. 
Of course, we could expect that the SPE possesses a loop soliton solution, since 
this equation is associated with a WKI-type spectral problem, like the equation 
studied in (jj] does, which describes the propagation of a loop along a stretched 
rope. The electric field in an optical fiber, however, is not a stretched rope, and 
the SPE solution for an ultra-short light pulse we are looking for must be neither 
multi- valued nor singular one. 

It is quite obvious that any solution of the SGE, which contains an asymptotically 
free kink (or antikink) at large t, corresponds via to a solution of the SPE 

containing a loop (resp., antiloop). Let us illustrate this with the two-kink solution 



z = —4 arctan 
and the kink-antikink solution 

z = —4 arctan 



m cosh ip 
n sinh (f) 

m sinh ip 
n cosh d> 



of the SGE © [H], where 

<f> = m(y + t), ip = n(y — t), n = \J m 2 — 1, 



(14) 
(15) 
(16) 



ANTON SAKOVICH AND SERGEI SAKOVICH 





Figure 2. The two-loop solution u(x,t) lO of the SPE |TJ| with 
to = 1.2: (A) t = -1.6; (B) t = 0.9. 




FIGURE 3. The loop-antiloop solution u(x,t) JTS1) of the SPE 
with to = 1.2: (A) t = -1.6; (B) t = 0.9. 



to is a real parameter, to > 1. Applying the transformation (JHJ) © to the ex- 
pressions Ijl4(l and Ijl5|l , and using the Mathematica system ^0] to do the required 
symbolic integration and simplification, we get, respectively, the two-loop solution 



Amn 



m cosh ip cosh (/) + n sinh ip sinh i 



x = y + 2mn 



to 2 cosh 2 tp + n 2 sinh 2 
to sinh 2-0 — n sinh 20 



(17) 



m 2 cosh 2 ip + n 2 sinh 2 
and the loop-antiloop solution 

to sinh tp sinh + n cosh -0 cosh 



u = 4mn- 



y + 2mn 



to 2 sinh 2 %p + n 2 cosh 2 q 
to sinh 2-0 — n sinh 20 



(18) 



to 2 sinh 2 ip + n 2 cosh 2 

of the SPE Q), with the notations i|16|) . Since these solutions (II 7|) and (|18fl . shown 
in Figures |2] and |21 are manifestly multi- valued, they cannot represent light pulses 
in optical fibers. 

It is easy to find out why and when the transformation (JHJ @: being applied 
to a solution of the SGE, yields a solution of the SPE with a singularity. From the 
relations JHJ and @ we have 



u x (x,t) = tan z(y, t). 



(19) 
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Then it immediately follows from Ijl9(l that a target solution u(x,t) of the SPE 
can be free from singularities only if the corresponding source solution z(y,t) of 
the SGE nowhere reaches any of the values (fc + i) it with k = 0, ±1, ±2, .... Of 
course, the solutions 1|12|) . I|14|) and (|15(l which we have already used, as well as 
all other solutions of the SGE which contain asymptotically free kinks or antikinks 
at large t, do not satisfy this requirement. For this reason, if we want to obtain 
nonsingular solutions of the SPE, we have to use those solutions of the SGE which 
correspond to bound states of kinks and antikinks. 

Let us take the breather solution of the SGE which is known to be the bound 
kink-antikink state 

/ msin?/> \ 

z = -4 arctan -— , (20) 



i cosh ( 



where 



<p = m(y + t) , ip = n(y — t), n = y 1 — to 2 , (21) 

to is a real parameter, < to < 1. Applying the transformation JHJ I© to the 
expression (|20[1 , we get the following solution of the SPE : 

to sin ip sinh <p + n cos tp cosh (f> 

u = Amn ;— J = ■ 

to 2 sin tp + n 2 cosh (j) 

m sin 2-0 — n sinh 2<p 
x — y + 2mn ;— = g — i 

to 2 sin if) + n 2 cosh 

with the same notations Q2ip. Note that we used the breather solution of the SGE 
in the form l|20() simplified by a Lorentz transformation and shifts of y and t, and 
the constant of integration in x (1221) has been fixed. If necessary, one can generalize 
the solution u(x, t) 1221) by the scale transformation and shifts of x and i. 

The solution l|22ll of the SPE is the central result of this paper, and hereinafter 
it will be referred to as the pulse solution. This name does not mean, however, that 
the expressions (|22|) represent a nonsingular pulse for any value of the parameter 
to, < to < 1. The function z(y, t) l|20() satisfies the condition — ^ < z < -| for all 
values of y and t only if to < m cr , where 

TO cr = sin f fa 0.383. (23) 

Therefore, when the value of to reaches or exceeds the critical value m CI (|23|) , the 
singularities u x — > ±oo appear in the function u(x, t) (|22|1 due to the relation (|19|l . 
Furthermore, in the overcritical case, when to > ui CY , the pulse solution (1221) is 
multi- valued, as shown in Figure 0] 

In the undercritical case, when to < to ct , the pulse solution 1221) represents a 
single- valued nonsingular pulse or, more precisely, a wave packet. Its typical shape 
is shown in Figure [SJ The envelope curve of this wave packet is determined by the 
hyperbolic functions in the expressions (1221) and moves from the right to the left with 
unit speed, whereas the oscillatory component is determined by the trigonometric 
functions in Q22JI and moves from the left to the right, also with unit speed. The 
smaller is the value of the parameter to in the pulse solution 11221) , the larger is the 
number of oscillations in the represented wave packet, as shown in Figure |H1 When 
the value of to is small, the pulse solution (|22ll can be approximated as 

u pa 4tocos(?/ — t) sech(m.(y + <)), x~y, (24) 
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FIGURE 4. The pulse solution u(x,t) (J2U of the SPE JTJ with 
to > to ci .: (A) m = 0.5, t = -0.6; (B) to = 0.8, t = -0.9. 





FIGURE 6. The pulse solution u(x 1 t) of the SPE with 
to < to ct , at t = 0: (A) to = 0.1; (B) m = 0.3. 



and the shape of the SPE pulse is similar to that of the NLSE soliton: the sech- 
shaped envelope modulates the cos-shaped wave (see Figure EJV) . On the other 
hand, very short pulses of the SPE, like the one shown in Figure|Hf3, correspond to 
near-critical values of to, for which the approximation (|24|l does not work fine and 
one has to use the exact expressions (122t . The ultra-short pulses represented by the 
pulse solution 111' 11 with near-critical values of to can be as short as approximately 
three cycles of their central frequency. 

In the critical case, when to = m cr = sin the pulse solution (|22|l still remains 
single- valued, but the singularities u x — > ±oo appear in it. This critical wave packet, 
or the shortest single- valued pulse of the SPE which we were able to find using the 
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FIGURE 7. The pulse solution u(x,t) of the SPE JTJ with 
m = m cr fa 0.383: (A) t = -0.7; (B) t = 0.5. 

transformation JSJ)-© between the SPE (Q) and the SGE is shown in FigureQ 
At present, we do not know whether the SPE can describe a shorter wave packet. 
It may well be that not, and that a modification of the model equation is necessary 
for extremely short optical pulses in nonlinear media, like the one proposed very 
recently in 

In conclusion, we can formulate the main result of the present paper as follows. 
We have found an exact nonsingular solitary wave solution of the Schafer- Wayne 
short pulse equation. This pulse solution was derived from the breather solution 
of the sine-Gordon equation by means of a transformation between these two inte- 
grable equations. The ultra-short wave packets represented by the obtained pulse 
solution can be as short as approximately three cycles of their central frequency. 
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